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Abstract 

We report on a study of the thermal dust emission of 
the circumstellar envelopes of a sample of Class sources. 
The physical structure (geometry, radial intensity profile, 
spatial temperature and spectral energy distribution) and 
properties (mass, size, bolometric luminosity (Lbol) an d 
temperature (Tboi), and age) of Class sources are derived 
here in an evolutionary context. This is done by combining 
SCUBA imaging at 450 and 850 /im of the thermal dust 
emission of envelopes of Class sources in the Perseus and 
Orion molecular cloud complexes with a model of the enve- 
lope, with the implementation of techniques like the black- 
body fitting and radiative transfer calculations of dusty 
envelopes, and with the Smith evolutionary model for pro- 
tostars. The modelling results obtained here confirm the 
validity of a simple spherical symmetric model envelope, 
and the assumptions about density and dust distributions 
following the standard envelope model. The spherically 
model reproduces reasonably well the observed SEDs and 
the radial profiles of the sources. The implications of the 
derived properties for protostellar evolution are illustrated 
by analysis of the I/boi, the Tboi, and the power-law index p 
of the density distribution for a sample of Class sources. 
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1. Introduction 

The earliest evolutionary phase of star formation, in which 
a deeply embedded protostar is known to exist, is the so- 
called Class stage (Andre et al. 1993). Since protostars 
in this phase (Class sources) are highly obscured by ex- 
tended dusty envelopes, they can be observed mainly in 
the far-infrared to millimetre wavelength range. Specif- 
ically, a Class source consists of a central protostellar 
object surrounded by an infalling dusty envelope and a 
flattened accretion disk. Many (presumably all) Class 
sources are associated with strong molecular flows, usually 
more energetic than those from Class 1 sources (Davis & 
Eisloffel 1995). 

Since Class sources will determine most of the early 
evolutionary history of the protostars as they evolve to- 
ward main-sequence stars, it is crucial to determine the 



physical conditions within the envelope of Class objects. 
However, because of detection difficulty, constraining phys- 
ical properties of Class sources is a challenging task. Nev- 
ertheless, continuum submillimetre observations of these 
sources let us to detect dust emission of the massive cir- 
cumstellar envelopes, and provide a powerful tool for con- 
straining the distribution of matter (Adams 1991). 

Independent methods of detection of Class sources, 
and techniques and useful ways to help to the interpreta- 
tion of properties of protostars have been developed and 
used during the last years (e.g. comparation of submil- 
limetre with near infrared surveys of a star forming region 
(Stanke et al. 2000), interpretation from the dust emission 
(Shirley et al. 2000), blackbody fitting (Chini et al. 2001), 
radiative transfer codes (Wolf et al. 1999), evolutionary 
models for protostars (Smith 2002)). In this work, we com- 
bine them and develop an approach in order to constrain 
the structure (geometry, radial intensity profile, spatial 
temperature, and spectral energy distributions (SEDs)), 
and physical properties (mass, size, bolometric luminosity 
(Lboi) and temperature (Ibd), and age) of Class sources. 
Specifically, we combine the information from the dust 
emission with a physical model of the extended envelope, 
the implementation of techniques like the blackbody fit- 
ting and the ID continuum self-consistent MC3D radiative 
transfer code (by Wolf et al. 1999), and the Smith evolu- 
tionary model for protostars (by Smith 2002). 

To describe and study stellar evolution, the Hertzsprung 
-Russell Diagram (HRD) has become the main way. How- 
ever, to study protostellar evolution, the Tboi-Lboi dia- 
gram (BLTD, bolometric luminosity-temperature diagram) 
was proposed as direct analogue to the HRD (Myers & 
Ladd 1993; Myers et al. 1998). There, Class sources are 
characterised by the lowest Tboi values of all protostars. 
Although the BLTD provides a useful way to study the 
evolution of protostars for all evolutionary phases of rele- 
vance for star formation (Myers et al. 1998), nevertheless, 
the Smith scheme predicts the evolution of Tboi and Lboi 
as an envelope of initial mass Mq and temperature T dis- 
sipates on a timescale t and forms a star. 

2. Theory 

How can we describe the protostellar emission from the 
envelope? In order to facilitate the interpretation of sub- 
millimetre observations, we adopt the standard envelope 
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model (Adams 1991) to keep the problem as simple as pos- 
sible. The circular symmetry of the observed emission (see 
Rcngcl 2004), and the lack of significant internal struc- 
ture justifies the simplicity of a spherical model case. If 
the emission is optically thin, the observed intensity for a 
spherically symmetric protostellar envelope at an impact 
parameter b is given by Eq. 1 assuming a single power-law 
index of the opacity [3. 
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I v {b) = 2K v B v [T d {r)]p{r) 
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ro is the outer radius, p the density, Td the dust tempera- 
ture, k„ the opacity of the dust grains, and B„[Td(r)] the 
Planck function at dust temperature Td- 

The current theory suggests that the density and tem- 
perature distributions are power-law indices (p and q, re- 
spectively, see e.g. Adams 1991). If the emission is in the 
Rayleigh-Jeans limit and if ro3> b, Eq. 1 can be approxi- 
mated to I v (b)/I V {Q) = (b/bo)~ m (where m is the power- 
law index of the observed intensity profile, and bo is the 
normalization factor to the peak emission). 

3. Observations 

To investigate the physical structure, processes, and prop- 
erties of Class sources, six star forming regions were 
observed in the Perseus and Orion molecular cloud com- 
plexes at 450 and 850 /im with the Submillimetre Common- 
User Bolometer Array (SCUBA) camera at the James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea, Hawaii. 
The target regions were L1448, L1455, NGC 1333, HH211, 
L1634, and L1641 N. A distance of 300 pc for L1448 is as- 
sumed according to Bachiller & Cernicharo 1986, 350 pc 
for L1455, NGC 1333 and HH211 (Chandler & Richer 2000, 
and Herbig & Jones 1983), 460 pc for L1634 (Bohigas et 
al. 1993), and 390 pc for L1641 (Anthony-Twarog 1982). 

As an example, the resulting combined SCUBA maps 
at 450 and 850 pm of L1455 overlaid with isophotal con- 
tour plots are shown in Fig.2] 

From these continuum SCUBA maps, 36 submillime- 
tre sources (> 3cr) are detected (26 in Perseus and 10 in 
Orion). Some of them are extended, and many contain 
multiple condensations, as well as extended diffuse fea- 
tures. Twelve objects are reported for the first time (six in 
Perseus and six in Orion) , and fifteen sources are analyzed 
here. Additionally, the first detailed submillimetre obser- 
vations of the two regions NGC 1333 South and L1641 N 
have been obtained and are discussed in Rengel (2004), 
unveiling objects previously unknown. The survey also re- 
veals a large number of small groups (Rengel et al. 2003). 

3.1. Derivation of parameters from the emission 

Physical parameters like the dust submillimetre spectral 
index 0450/850 (this value is related directly to the dust 
emissivity exponent (3 by (3~a - 2), gas and dust masses, 
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Figure 1. Maps of the region L14-55, overlaid with con- 
tours. The top map is taken at 850 pm and the bottom 
map at 450 pm. In L 1455 the protostars RNO 15 FIR and 
RNO 15, and the pre-stellar cores SMS I and II are ob- 
served. Contour levels are in log scale with step size 1.15 
from 0.0005 Jy/beam at 850 pm, and step size 1.25 from 
0.0035 Jy/beam at 450 pm. The filled circles at the bottom 
left indicate the main beam size. 



and sizes of the envelopes are derived. The dust submil- 
limetre spectral index cxx 1 /\ 2 is defined by 



(2) 



log(X 2 /Xi) 

The mean value of 0450/850 f° r these objects is 2.8 ± 
0.4, implying that the regions investigated here are quite 
cold (~10 K). Such low dust temperatures have for exam- 
ple also been found for isolated globules in the Perseus 
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region (see e.g. Myers & Benson 1983), suggesting an ex- 
ponent ~ 1 for the dust grain opacity, K v . The failure 
of the Rayleigh-Jeans approximation could increase the 
value of (3. 

An aspect ratio was defined for each of the sources by 
dividing the semi-major and minor axis lengths as mea- 
sured typically at the 97% contour. We find that the ob- 
served aspect ratios of the protostellar emission maps lie 
in the range 1.0-2.0 with an average value of 1.3 ± 0.1, 
suggesting that the envelopes of Class sources can be 
described approximately as spherically symmetric. Never- 
theless, a departure of this symmetry, i.e., from an aspect 
ratio of 1 , in some cases is found. This mismatch could be 
produced by several factors like magnetic fields of energy 
densities sufficiently large to influence the object structure 
and cause flattening along field lines, the bipolar outflow, 
rapidly rotating structures, or high levels of angular mo- 
mentum (Rengel 2004). 

The average value of the gas and dust mass of the sam- 
ple, based on the thermal emission from the dust according 
to Hildebrand (1983), is 2.5 ± 0.6 M Q (typical aperture 
45"). In order to derive physical properties of the sources 
from the observed azimuthally averaged radial intensity 
profiles of the thermal emission from the dust, like the 
power-law index of the temperature distribution as a func- 
tion of the radius T(r) cx r~ q , q, p, and the envelope sizes, 
the standard envelope model (Adams 1991) was adopted. 
Assuming that the observed sources can be fitted by a 
power-law intensity distribution I v (b)/Io = (b/bo)~ m as 
Class0/1 sources, the average value of <m> is 1.74±0.02 
at 850 /im and 1.50 ± 0.02 at 450 ^m. Following the an- 
alytic estimate for the intensity (Adams 1991), these val- 
ues of <m> correspond to the expected ones for very 
young objects: theoretical models as well as numerical 
simulations of the collapse of an isothermal sphere (e.g., 
Shu 1977; Larson 1969; Penston 1969) also predict power- 
law indices of the density distribution, p(r) oc r~ p , with 
p < 2.0. Considering m—p + q -1, the first two order terms 
of the equivalent power-law index \x of the observed inten- 
sity (Adams 1991), the values q= 0.42 ± 0.04 as well as 
p= 2.1 ± 0.1 and 2.3±0.1 (at 450 and 850 fim, respec- 
tively) are found. In Fig.|2 the normalized radial pro- 
files [I v (b)/I u (0)} (being I„(Q) the observed peak flux) of 
RNO 15 FIR as an example are plotted versus b. 

To estimate the source size, an outer diameter d corre- 
sponding to the measured Half Power Beam Width (HPBW) 
size is adopted. The observed sources are surrounded by 
extended envelopes, having typical sizes of 1500-6000 AU 
(at 450 fim) and 4000-9000 AU (at 850 fim). 

4. Modelling 

4.1. Blackbody fitting 

Accurate SEDs are computed for nine sources combin- 
ing existing multi-wavelength surveys with the submil- 
limctre data obtained in Sect. 3 (Rengel 2004), and fitting 
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Figure 2. Radial profiles of RNO 15 FIR at 850 fim (top 
panel) and 450 fim (bottom panel). The normalized inten- 
sity is plotted as a function of impact parameter, b(AU). 
The bold solid lines represent the fitted slope m. The hori- 
zontal dashed lines represent the 1 a noise level at the edge 
of the map. The beam profiles are shown as dashed lines. 
Derived diameters d of the source (in AU) are given in the 
lower parts. 

a modified blackbody curve to the fluxes (Froebrich et 
al. 2003). This method computes a grid of different black- 
body curves of the form J2J in which the temperature, the 
optical depth at 100 fim and f3 are varied (15 K <T< 
80 K; 0.0 < [3 < 3.0; 0.09 < t wo < 11.4). Specifically, 
curves of the form Q are convolved with the filter curves 
of the used filter bands (Eq.EJ), and compared to the data- 
points. The model that best fits the observational data is 
determined. We derived the bolometric temperature and 
luminosity, size of the envelope, and sub-millimetre slope 
for each source. 

In local thermodynamic equilibrium, the solution of 
the equation of transfer for an isothermal object that sub- 
tends an effective solid angle EJ7 is given by the emergent 
flux density S: 

S[Jy] — Sfi • (l — e~ r ) ■ B(\, T d ). (3) 

B(X,T) denotes the Planck function for a dust tempera- 
ture Td, and r the dust optical depth. The optical depth 
can be expressed by 

/ A \ _/3 

T = Tl0 ° ' Tori ' W 

\100^my 
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A is in /Ltm, the optical depth at 100 /mi (tiqo) is a free 
parameter, and (3 is the submillimctre slope of the SED 
(/3=0 corresponds to a blackbody). 

But we do not observe the emergent flux density S\ 
directly, we observe Sf, the emergent flux detected across 
a filter /, which has a filter transmission curve 7/ (A). Sf 
is calculated by 

£0 J (1 - e- T ) B(X,T)T f (X)dX 

Sf [Jy] = — ^ 5S (5) 

/ T f (X)dX 

for each filter / separately using the filter transmission 
curves Tf(X). For details about how the calculations of 
the maximum likelihood L were performed, see Froebrich 
et al. (2003). 

The SEDs revealed that the sample consists of cool 
objects (Tboi ranges of ^27-50 K) which have Lboi of ~4- 
85 L/0. Thus, by means of the sub- millimetre to bolometric 
luminosity ratio, L1448NW, L1448C, RNO 15 FIR, NGC 
1333 IRAS 1, NGC 1333 IRAS 2, HH211-MM, L 1634, L 16 
41 N, and L1641 SMS III are Class sources. For detailed 
computations of the associated errors to each determined 
parameters, see Rengel (2004). In Fig.0] (bottom, short 
dashed lines), best obtained fit to the data-points using 
our approach and overplotted the combined data of RNO 
15 FIR are shown as an example. 

Because of the IRAS 12 and 25 /an points are usually 
far above the fit to the SED (e.g. Chini et al. 2001; Barsony 
et al. 1998), and of the unclear quantification of a higher 
bolometric temperature value caused by filter leaks, out- 
flow/dust interactions or envelope emission if these IRAS 
points are included on SED computations (^4 K Froe- 
brich et al. 2003), these points are excluded on the target 
list and on the blackbody fits performed here. In Fig.0] 
the blackbody curve and the data-points show significant 
deviations at the Infrared Astronomy Satellite (IRAS) 12 
and 25 ^m points for the Class sample. If these data- 
points had correct flux measurements, the sources would 
exhibit an excess of mid-infrared emission which could be 
caused by ongoing outflow/dust interactions. Because it is 
likely that these IRAS points are incorrect in some cases, 
further mid-infrared observations of these objects are nec- 
essary in order to investigate the nature of this possible 
excess. 

4.2. MC3D RADIATIVE TRANSFER CODE 

Results from the standard theory provide the first di- 
rect insights into observable estimations (e.g. p and q, 
Sect. 3.1). Nevertheless, the temperature profile will di- 
verge from a single power-law q as the envelope becomes 
optically thick at the primary wavelengths of energy trans- 
port (inner portion of the envelope) (e.g. Shirley et al. 2000) 
Here it becomes necessary to calculate the temperature 



distribution self-consistently by implementation of a ra- 
diative transfer code. 

Detailed physical interpretation of the internal struc- 
ture of nine objects is realized by intensive computer ra- 
diative transfer modelling. In particular, dust temperature 
distributions, SEDs, and intensity maps are derived self- 
consistently by the MC3D code (Wolf et al. 1999). The 
basic model assumed is: a central heating source that is 
embedded within a spherical symmetric dust envelope of 
mass M cnv . The central source is characterized by a stel- 
lar luminosity (£*), and an effective temperature (T*). 
The dust envelope is specified by an outer radius (i? ut), 
a density power-law index p, an inner radius of the enve- 
lope i?i n , which is set by the dust destruction radius i? su b 
(also called sublimation radius) adopted to a temperature 
of 1500 K, the dust destruction temperature of graphite 
(e.g. Adams & Shu 1985, which occurs at a temperature 
of 1500 K; Myers et al. 1998). The central source lumi- 
nosity is the sole source of luminosity for the envelope, 
so L„=Lh \. Fig.|3]shows a scheme of the adopted model. 
The validity of this model, and physical parameters of the 
sources (e.g. envelope masses, density distributions, sizes, 
and sublimation radius) are derived by finding the consis- 
tency between observed and modelled SEDs, and radial 
profiles (Sect 4.3; Rengel et al. 2004). 

The dominant opacity source in the radiative trans- 
fer models is dust grains. The characteristics of grains 
in Class sources are unknown. However, graphite and 
silicate grains have been proposed as likely consistent of 
Galactic dust. The grain distribution adopted here con- 
sists of small interstellar grain size 0.005-0.25 /im, with 
mixing of astronomical silicate and graphite (relative abun- 
dances: 62.5 % & 37.5 % (respectively)), and the Draine 
& Lee (1984) optical constants to describe grain proper- 
ties. The adopted model assumes that the dust properties 
are not a function of radius in the envelope. The gas-to- 
dust mass ratio is adopted to be 100:1, whereby we use a 
dust grain density of 3.6 g cm~ 3 . The wavelength grid (60 
wavelengths) are chosen to cover the relevant range. 

In order to derive the physical structure of the enve- 
lope, three steps are carried out: the first step constrains 




Figure 3. Basic model assumed in the MC3D code. 
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the spatial temperature structure of the circumstellar en- 
velope, the second step determines the simulated SED, 
and the third step, derives the emission, which provides 
the intensity structure. These last observables lead to test 
the model. At this point it is necessary to keep in mind 
that the values of the inferred parameters are valid only 
within the framework of the model assumptions consid- 
ered here. In Fig.0] the temperature profile, and best sim- 
ulated and observed SEDs for RNO 15 FIR are shown as 
an example. 

The central source should radiate from a specific radius 
r, where the layer changes between the optically thick and 
the thin regime. This radius r is called "photosphcric ra- 
dius" or "effective photosphere", which depends on the 
wavelength and the molecule species. Crudely, the ob- 
served emission arises from a "characteristic" layer at t\=1 
(Myers et al. 1998) or 2/3 (Kenyon et al. 1993; Hart- 
mann 1998). The temperature distribution changes its char- 
acter inside the photospheric radius, where the envelope is 
optically thick; trapping of radiation causes the tempera- 
ture gradient to become steeper than in the optically-thin 
case. To calculate the photospheric radius, r=l and the 
spherically symmetric optically thick envelope of Myers et 
al. (1998), with an inner radius R ln , a density p ln and an 
outer radius i? ut , are adopted here for simplicity. Follow- 
ing this model, the photospheric radius r is given by 

r T= i S3 -Ri n Tl2, (6) 

T12 being the optical depth from i?j n to i? ou t- 

Following the Myers et al. (1998) model, a sublima- 
tion radius of 3-5 AU at 450 /im, r=l and p=2 (~10 AU 
for p=1.5) is estimated, which is roughly the same as the 
radius of the photosphere. 

As result of the modelling, it is found that the physi- 
cal structure of a Class envelope is characterized by the 
gas and dust temperature distribution T(r) oc r~~ q , and 
the density distribution p(r) oc r~ p , with q=0A and p 
in the range of 1.5-2. In the inner part of the envelope 
(10 AU), T(r) departs from the single-power law index of 
0.4, resulting in a steeper inner temperature profile. This 
phenomenon can be produced probably by thermal con- 
vection in the inner envelope and causing zones of dust 
destruction by evaporation when the envelope becomes 
optically thick (e.g. Lin & Papaloizou 1980). T accounts 
for the effects of this heating. 
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Figure 4- Temperature profiles for RNO 15 FIR performed 
using the MC3D radiative transfer code (Wolf et al. 2003) 
(top). Best simulated and observed SEDs (bottom, short 
dashed and dotted lines, respectively). The filled squares 
represent the measured 450 and 850 [im fluxes from the 
SCUBA observations presented in Sect. 3. The IRAS 12 
and 25 fim points are excluded from the fit of the SED. De- 
tails regarding the summarized target list at several wave- 
lengths from the literature are given in Rengel (2004). For 
the model results, see Table 1. 



4.3. Comparison of models with observations 

We estimate the parameters for the Class sources by 
comparing the observed SEDs obtained by blackbody fit- 
ting fSect. l4~T|l and the simulated SEDs from the envelope 
fitting procedure fSect. l4~2*)l . Following the same technique 
as Kenyon et al. 1993, the best model for a given observa- 
tional SED is determined by minimizing the normalized 
sum of the weighted squared residuals, 



S 2 = 



IVii^Foji Fm,i) i Fq : : 



where F D ^ is the observed flux at wavelength i, F m ^ the 
scaled flux of the model, and u>i a weighting factor 1 . Given 
the uncertainties in estimating observational errors and 
the limitations of the models, a more complicated fit- 



1 a weighting factor is empirically given to the consistency 
between the modeled and observed SED, with dependence on 
the wavelength. Reasonable weights are given from the uncer- 
tainties in each region of the SED. 
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ting treatment is unjustified. Because of the uncertain- 
ties in opacity, the difficulties in matching observational 
beam size with model results, and the possible contribu- 
tion of disks to the submillimctre region (c.f. Beckwith et 
al. 1990; Butnet et al. 1991; Keene & Masson 1990), we 
used weights Wi=l for wavelengths 60 fim < A < 300 /xm, 
and chose u>i=0.5 for 300 /mi < A < 1000 /im. Data-points 
taken at wavelengths shorter than 40 /Ltm are excluded in 
the comparation calculations because they are expected 
to be optically thick and therefore very sensitive to asym- 
metric geometries (e.g. outflow cavities and flattened en- 
velopes), which are currently unable to be modelled. Be- 
cause of the spectral range considered here does not extend 
much over long wavelengths, data-points taken at wave- 
lengths longer than 1000 fim are rejected. 

As it is claimed previously, a comparison between best 
simulated and observed SEDs, and between simulated and 
observed radial profiles led us to check the reliability of 
the model. In order to create the simulated radial profile, 
the envelope emission at 850 fim is calculated as an image 
from the MC3D code for each source, and convolved by the 
SCUBA beam. Fig.[S] shows the observed and simulated 
radial profile of RNO 15 FIR as an example. 

We conclude from our comparison that the assumption 
of spherical symmetry in the envelope is a reasonable ap- 
proximation to study the envelope of Class sources. Nev- 
ertheless, deviations from this geometry between observed 
and modeled SEDs at short wavelengths, and between ob- 
served and simulated radial profiles of relatively spheri- 
cal Class sources become extremely prominent (e.g. for 
RNO 15 FIR and NGC 1333 IRAS 1) (see Rengel 2004). In 
addition, small differences between observed and simu- 
lated radial profiles are observed for HH 211-MM and RNO 
15 FIR (aspect ratios of 1.2 and 1.3 at 850 /im, respec- 
tively). This discrepancy could be due to the exclusion of 
a non-spherical geometry, of cavities produced by molec- 
ular outflows, or of a flattened disk in the model; it could 
also be due to specific assumptions considered here, e.g. a 
single point-like nature of the sources (e.g. it is uncertain 
that HH211-MM and RNO 15 FIR are point-like sources 
(Eisloffel et al. 2003; Davis et al. 1997, respectively)), and 
a constant dust opacity with radius in the model. 

4.4. BOLOMETRIC LUMINOSITY-TEMPERATURE DIAGRAM 

The observed properties of Class sources can be under- 
stood more thoroughly if it is assumed that the objects 
evolve in time. Using the Smith protostellar evolutionary 
scheme (Smith 2002), and constructing the Izboi-^boi di- 
agram for the sample, the age, the infall rate, and the 
envelope mass as a function of time are estimated for the 
detected Class sources (Rengel 2004). In particular, com- 
paring the positions of the objects with an evolutionary 
model age and mass of the envelope of each source are 
derived (Fig.|Hl filled triangles). 



The tracks in Fig. are derived by combining the Uni- 
fication Scheme, as reviewed by Smith (2000, 2002), with 
the framework for protostellar envelopes presented by My- 
ers et al. (1998). 

For the sample, the resulting model-dependent ages of 
Class sources are in a range from 1 to 3 x 10 4 yr. The 
evolution of these sources is traced by Tboi , which typically 
increases systematically between 25-50 K. The resulting 
bolometric luminosities span a range between 5-85 L Q . 

Following the Smith scheme (with an accretion rate 
not constant in time, Smith 2002), the accretion rate M acc 
at the earliest Class stage increases exponentially for a 
short interval of time (with an accretion peak in 1.5xl0 4 
yr), and then declines later on (Rengel 2004). 

Is there a correlation of the power-law index p with 
time? We plot the estimated values of p as a function 
of t. Fig. suggests a density structure of p(r) <x r~ 2 at 
younger ages, evolving to p(r) cx r -3 / 2 at later times. In 
order to investigate this with a larger spread, a Class 1 
sample is mandatory. 

5. Results 

Table 1 lists the characteristics of the best-fitting model 
parameters for the sample of sources analyzed in detail 
here. 

The main results of this study are: 1) From the con- 
tinuum SCUBA maps, 36 submillimetre sources (> 3cr) 
are detected. Some of them are extended, and many con- 
tain multiple condensations, as well as extended diffuse 
features. Twelve of these objects are reported here for 
the first time, 15 sources are analyzed, and nine Class 
sources are studied in detail. 2) An average value of the 
aspect ratio of the sample of 1.3 ± 0.1 is found, suggest- 
ing that the envelopes of Class sources can be described 
approximately as spherically symmetric. Nevertheless, a 
departure of this symmetry in some cases is found. This 
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Figure 5. Simulated and observed radial profiles as a func- 
tion of the impact parameter b(AU) of RNO 15 FIR at 850 
fim (short dashed and solid lines, respectively) . The dotted 
lines represent the 1 a noise level at the edge of the map. 
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Figure 6. The BLTD for the objects analyzed here (filled 
triangles), the Class data from Andre et al. (2000) (open 
squares), from Froebrich et al (2003) (thick Xs), and from 
Stanke 2000 (open diamonds). Protostars evolve from 
right to the left. Three tracks for the final masses of 0.2, 
1 and 5 M Q are displayed. The vertical dotted lines on the 
tracks mark the model ages of 20, 30, 40, and 50 thou- 
sand years (tracks derived by combining the Unification 
Scheme (Smith 2000; 2002) with the framework for proto- 
stellar envelopes presented by Myers et al. (1998)). 



mismatch could be produced by several factors like mag- 
netic fields of energy densities sufficiently large to influence 
the object structure, the bipolar outflow, rapidly rotating 
structures, or high levels of angular momentum. 3) The 
sample consists of cool objects (Tboi ranges of —27-50 K) 
which have Lboi of —4-85 L Q . 4) The physical structure 
of a Class envelope is characterized by the gas and dust 
temperature distribution T(r) oc r~ q , and the density dis- 
tribution p(r) oc r _p , with q=0A and p in the range of 
1.5-2 (p—2 at younger ages and p—1.5 at later times). In 
the inner part of the envelope (10 AU), T(r) departs from 
the single-power law index of 0.4, a phenomenon caused 
probably by thermal convection in the inner envelope. 5) 
From the thermal emission of the dust and following an 
emissivity model (Hildebrand 1983), the average value of 
the gas and dust mass of a sample based on 15 submil- 
limetre selected objects is 2.5 ± 0.6 M©. Examination of 
the radial profiles of the sample shows that the objects 
are surrounded by extended envelopes. Typically, the sizes 
are -1500-6000 AU (at 450 /zm) and -4000-9000 AU (at 
850 /um). 6) A sublimation radius of 3-5 AU at 450 /im, 
t=1 and p~2 (—10 AU for p=1.5) is estimated, which is 
roughly the same as the radius of the photosphere (be- 
cause the gas that is accreted close to the central object 
is optically thick, the central source should radiate from 
this specific radius, where the layer changes between the 
optically thick and the thin regime). 7) The blackbody 



_j i I i i i i I i i i i I i i i 

1.5xl0 4 2xl0 4 2.5xl0 4 
age [yr] 

Figure 7. The power-law index p of the sample as function 
of age listed in Table 1. Detailed computations show that 
the errors of p are less or equal to 0.1. 

curve and the data-points show significant deviations at 
the Infrared Astronomy Satellite (IRAS) 12 /im and 25 ^m 
points for the Class sample. If these data-points had cor- 
rect flux measurements, the sources would exhibit an ex- 
cess of mid-infrared emission which could be caused by 
ongoing outflow/dust interactions. Because it is likely that 
these IRAS points are incorrect in some cases, further 
mid-infrared observations of these objects are necessary 
to investigate the nature of this possible excess. 8) The 
assumption of spherical symmetry in the envelope model 
is reasonable. Nevertheless, deviations from this geometry 
determined by comparations between observed and mod- 
eled SEDs at short wavelengths, and between observed 
and simulated radial profiles of relatively spherical Class 
sources become extremely prominent. This discrepancy 
could be due to the exclusion of a non-spherical geom- 
etry, of cavities produced by molecular outflows, or of a 
flattened disk in the model; it could also be due to spe- 
cific assumptions considered here, e.g. a single point-like 
nature of the sources, and a constant dust opacity with 
radius in the model. 9) For the sample of Class sources, 
the model-dependent ages are in a range of 1-3 xlO 4 yr. 

6. Conclusions 

Detailed physical interpretation of the internal structure 
of nine observed Class sources was carried out by inten- 
sive computer modelling. A simple spherically symmetric 
model envelope, and assumptions about density and dust 
distributions following the standard envelope model repro- 
duced reasonably well the observed SEDs and the radial 
profiles of the sources. T(r) oc r~ 0A is a good approxima- 
tion for the sample. The radial temperature distribution as 
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function of distance, however, departed significantly from 
the optically-thin assumption, an observational derivation 
of a single-power law q of 0.4 for radii < 10 AU. Mod- 
elling results indicated a density profile well described by 
a power-law between p=1.5 and 2, which is expected by 
all of the collapse models and numerical studies. A density 
structure of p(r) oc r~ 2 at younger ages, evolving to p(r) 
oc r~ 3 / 2 at later times was found. 

SED and radial profile fits constrained physical param- 
eters of the sources such as envelope masses, density dis- 
tributions, sizes and sublimation radii. Nevertheless, mid- 
infrared observations, and the inclusion of outflow and 
disk, and other geometries in the model are desirable to 
investigate in order to decrease the differences between 
the model and the data, or prove that data and model are 
different under several configurations. 
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Table 1. Best envelope fit results for the embedded sources. The temperature of the central object is set to 3500 K. 
X 2 quantifies the agreement between model and data. The age is estimated according to the model of Smith (2000) and 
given in l(r yrs. 



Object 


T bol [K] 


L* [Lq] 


Rsu b [AU] 


P 


Rout [AU] a 


,ge [x 10 3 yr] 


M cnv [M Q ] 


x 2 


L1448NW 


27 


6 


5 


2.0 


4000 


13.8 


2.8 


6.4 


L 1448 C 


40 


11 


5 


1.6 


3000 


19.5 


1.7 


2.6 


RNO 15 FIR 


45 


8 


3 


1.6 


3500 


20.7 


1.0 


1.8 


NGC 1333 IRAS 1 


52 


13 


3 


1.5 


4500 


21.9 


1.5 


8.6 


NGC 1333 IRAS 2 


48 


68 


3 


1.5 


7000 


26.1 


3.0 


0.1 


HH211-MM 


30 


5 


3 


1.9 


6000 


13.8 


4.1 


4.5 


L1634 


42 


19 


3 


1.6 


6000 


21.9 


2.8 


1.1 


L1641N 


45 


85 


3 


1.5 


8500 


25.2 


7.0 


1.0 


L 1641 SMS III 


50 


80 


3 


1.5 


10000 


27.6 


6.3 


0.6 



